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Abstract 

The behavior displayed by a quantum system when it is perturbed by a series of von Neumann measurements 
along time is analyzed. Because of the similarity between this general process with giving a deck of playing 
cards a shuffle, here it is referred to as quantum shujflmg, showing that the quantum Zeno and anti-Zeno 
effects emerge naturally as two time limits. Within this framework, a connection between the gradual 
transition from anti-Zeno to Zeno behavior and the appearance of an underlying Markovian dynamics is 
found. Accordingly, although a priori it might result counterintuitive, the quantum Zeno effect corresponds 
to a dynamical regime where any trace of knowledge on how the unperturbed system should evolve initially 
is wiped out (very rapid shuffling). This would explain why the system apparently does not evolve or decay 
for a relatively long time, although it eventually undergoes an exponential decay. By means of a simple 
■ working model, conditions characterizing the shuffling dynamics have been determined, which can be of 
help to understand and to devise quantum control mechanisms in a number of processes from the atomic, 
molecular and optical physics. 

Key words: Quantum shuffling. Quantum Zeno effect, Anti-Zeno effect. Measurement theory, 
von Neumann measurement, Markov chain 



1. Introduction 

The legacy of Zeno of Elea becomes very apparent through calculus, the pillar of physics. It is not difficult 
to find manifestations of his famous paradoxes throughout the subtleties of any of our physical theories [l[ . 
In quantum mechanics, for example, Zeno's paradox of the arrow is of particular interest, for it has given 
rise to what is now known as quantum Zeno effect (QZE), which constitutes an active field of research Q. 
As conjectured by Misra and Sudarshan Q, this effect essentially consists of inhibiting the evolution of an 
unstable quantum system by a succession of shortly-spaced measurements — a classical analog of this effect 
is the watched-pot paradox: a watched pot never boils U — , although, generally speaking, it could also be 
a system acted by some environment or even the bare evolution of the system if the latter is not described 
by a stationary state. The inhibition of the evolution of a quantum system, though, was already noted by 
von Neumann Q and others — an excellent account on the historical perspective of the QZE can be found 
in 0. From an experimental viewpoint, this effect was formerly detected by Itano et al. [y] considering the 
oscillations of a two-level system, a modification suggested by Cook 0] of the original theoretical proposal. 
Nonetheless, the first experimental evidences with unstable systems, as originally considered by Misra and 
Sudarshan, were observed later on by Raizen's group @, Indeed, in the second experiment reported by 
this group in this regard Q, it was also shown the possibility to enhance the system decay by considering 
measurements more spaced in time. This is the so-called (quantum) anti-Zeno effect (AZE) |1CH12|. 



In the literature, it is common to introduce the QZE and the AZE as antagonist, competing effects. In 
this work, however, we study their manifestation within a unifying framework, where they constitute the 
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two limiting cases of a more general process that we shall refer to as quantum shuffling. To understand 
this concept, consider a series of von Neumann measurements is performed on a quantum system, i.e., 
measurements such that the outcome has a strong correlation with the measured quantity, thus implying 
a high degree of certainty on the post-measured system state. This type of measurements provoke the 
system to collapse into any of the pointer states of the measuring device, breaking the time coherence 
that characterizes its unitary time-evolution. In other words, the continuity in time between the pre and 
post-measurement states of the system is irreversibly lost. This loss takes place even when the distance 
(in time) between the two states is so close that, in modulus, they look pretty much the same, due to the 
corresponding loss of the phase accumulated with time (which is a signature of the unitary time-evolution 
and therefore of the possibility to revert the process in time). Thus, consider the system is not stationary 
(regardless of the nature of the source that leads to such non-stationarity), decaying monotonically in time 
at a certain rate. A series of measurements on this system will act similarly to giving a deck of ordered 
cards a shuffle — hence the name of quantum shuffling — , affecting directly its coherence and modifying its 
natural (unperturbed) decay time-scale. As it is shown here, depending on the relative ratio between the 
natural time-scale and the shufflingly-modified one, the system decay can be either delayed or enhanced. If 
the shuffling frequency (i.e., the amount of measurements per time unit) is relatively low with respect to the 
system natural decay rate, the pre and post-measurement system states will be very different. This turns into 
a very fast decay due to the important lack of correlation between both states. Within the standard Zeno 
scenario, this enhancement of the decay corresponds to AZE. On the contrary, if the shuffling is relatively 
fast, the pre and post-measurement states will be rather similar, for the system did not have time enough 
to evolve importantly. The decay is then slower, giving rise to QZE. Now, as it is also shown, this inhibition 
of the system decay is only apparent: the fast shuffling gives rise to an overall exponential decay law at long 
times that makes the QZE to be sensitive to the total time along which the system is monitored. Thus, in 
the long term, one just finds out that the system evolution displays features typical of Markovian processes 
such as exponential decays (with relatively long characteristic times) and time correlation functions 
with the form of a Markov chain [14|. As a consequence, if the natural system relaxation goes as a decreasing 
power series with time (e.g., in systems with a regular system dynamics), in the QZE regime it is found 
that the perturbed system undergoes decays that fall below the natural decay after some time due to the 
exponential decay induced by the short-spaced measurement process. This unexpected behavior is usually 
missed and therefore unexplored, since in QZE scenarios it is more common to consider the mathematical 
limit rather than the physical one. 

In order to demonstrate the assertions mentioned above as clearly as possible, here we have considered as 
a working model the free evolution of a Gaussian wave packet. When it is not perturbed by any measurement, 
this non-stationary system undergoes a natural decay. That is, this decay is not bound to effects linked to 
the action of external potentials or surrounding environments, but only to the bare wave-packet spreading 
as time proceeds. From a time-independent perspective, this spreading is explained by the continuum of 
frequencies or energies (plane waves) that contribute coherently to the wave packet, which give rise to a 
non-stationary evolution with time; from a time-dependent view, it is just a diffraction effect associated with 
the initial localization (spatial finiteness) of the wave packet. In either case, this property together with 
the analyticity and ubiquity of the model (it is a prototypical wave function describing the initial state of 
atomic, molecular and optical systems) make of the free Gaussian wave packet an ideal candidate to explore 
the Zeno dynamics. As it is shown, specific conditions for the occurrence of both QZE and AZE are thus 
obtained in relation to the two mechanisms involved in its dynamics: its translational motion and its intrinsic 
spreading, which have been shown to rule the dynamics of quantum phenomena, such as interference [15] or 
tunneling [16]. More specifically, in order to detect QZE and AZE, the overlapping of the wave function at 
two different times must be non-vanishing. In this regard, therefore, if translation dominates the evolution, 
the correlation function will decay relatively fast and none of them will be observable. The analytical results 
here obtained, properly adapted to other contexts, may provide the physical insight necessary to understand 
more complex processes described by the presence of external interaction potentials or coupled environments 



17|,[18j. It is also worth stressing that, to some extent, the Zeno regimes found keep a certain closeness with 
the three-time domain scenario considered by Chiu, Sudarshan and Misra [l^l • 

The work is organized as follows. The dynamics of a free wave packet is introduced in Section [2l in 
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particular, the (analytical) behavior of its associated time-dependent correlation function. This will provide 
us with the basic elements to later on establish the conditions leading to QZE or AZE once the shuffling 
process induced by the succession of measurements will be introduced. In particular, the quantum shuffling 
effect is analyzed in Section [3] assuming the wave packet is acted by a series of von Neumann measurements. 
These measurements will be assumed to occur at equally spaced intervals of time and their action on the 
system will be such that the post-measurement state will always be equal to the initial one. This could 
be the case, for example, when considering projections (diffractions) through identical slits Finally, in 
Section!?] we summarize the main conclusions extracted from this work. 



2. Dynamics of a free wave packet 

2.1. Characteristic time scales 

Consider the initial state of a quantum system is described in configuration space by the Gaussian wave 
packet 

\]/q(x) = Aoe~'^^~^°''^^'^'^°'^*^"'^^^^°^^^', (1) 

where Aq = (27rtTg)~^/^ is the normalization constant, xq and po are, respectively, the position and transla- 
tional (or propagation) momentum of its centroid, and ctq is its initial spatial spreading. The time-evolution 
of this wave function in free space {V{x) = 0) is given [l5| by 

\I>j(a;) — j^^^-(x-xt)^ / Aaaat+ipa{x~xt)/h+iEot/h ^2) 

Here, At = (27r(T|)~^/^ is the time-dependent normalization factor; xt = X{) -\- v^t is the time-dependent 
position of the wave packet centroid, with vq — po/m being its speed; Eq — p'^/2m is the average translational 
energy, responsible for the time-dependent phase developed by the wave packet as time proceeds; and 
= Co [l -l- {iht/2ma'^)\ , with at = \£t \ = croy^l + {ht/2maQ)'^ being the time-dependent spreading of the 
wave packet. This spreading arises [l5[ from a type of internal or intrinsic kinetic energy, which can be 
somehow quantified in terms of the so-called spreading momentum, ps = h/2ao, an indicator of how fast 
the wave packet will spread out. This additional kinetic contribution becomes apparent when analyzing the 
expectation value of the energy or average energy for the wave packet, 

as well as in the variance. 




AE^J m - = J^J^ + ^. (4) 



(These two quantities are time- independent because of the commutation between Hamiltonian operator, H , 
and the time-evolution operator, U = e'*^^''.) From a dynamical point of view, the implications of this term 
in are better understood through the real phase of 

Six, t) = poix - xt) + (x - xtf + E,t^\ (tan)-i . (5) 

Putting aside the third and fourth terms in these expressions — two space-independent phases related to 
the propagation and normalization in time, respectively — , we observe that the first term is a classical-like 
phase associated with the propagation itself of the wave packet, while the second one is a purely quantum- 
mechanical phase associated with its spreading motion. Correspondingly, each one of these two motions 
leads to the two energy contributions that we find in ([3]). 

By inspecting the functional dependence of at on time, a characteristic time scale can be defined, namely 
T = 2ma^/h. This time scale is associated with the relative spreading of the wave packet, allowing us to 



distinguish three dynamical regimes in its evolution depending on the ratio between t and t 21|: 
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(i) The very-short-time or Ehrenfest-Huygens regime, t <^ r, where the wave packet remains almost 
spreadless: <7t ~ (To- 

(ii) The short-time or Fresnel regime, t t, where the spreading increases nearly quadratically with time; 
at~ao + (?iV8mVg)t2. 

(iii) The long-time or Fraunhofer regime, t ^ t, where the Gaussian wave packet spreads linearly with 
time: at « {h/2m(7o)t. 

By means of r we can thus characterize the dynamics of the Gaussian wave packet, although similar time 
scales could also be found in the case of more general wave packets provided that we have at hand their 
time-dependent trend — this is in correspondence with the time-domains determined by Chiu, Sudarshan 
and Misra for unstable systems [l^l- Keeping this in mind, consider the probability density associated with 
@, 

|vI/,(x)|2 = ^i=e-(---')'/2'^'. (6) 

Case (i) is not interesting, because it essentially implies no evolution in time. So, let us focus directly on 
case (ii), for which ^ reads as 
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TTO-n 



-{x-Xt) /2CTg 



(7) 



where the time-dependent factor in the argument of the exponential can be neglected without loss of general- 
ity (the exponential of such an argument is nearly one) . According to (O , the initial falloff of the probability 
density is parabolic and therefore susceptible to display QZE if a series of measurement is carried out at 



regular intervals of time [22, l23| provided that these time intervals are, at least. At < t (later on, in Sec- 
tion |31 another characteristic time scale, namely the Zeno time, will also be introduced). For longer time 
scales (case (iii)). 



Accordingly, for distances such that the ratio {x — Xt)/t remains constant with time (remember that the 
spreading is now linear with time), the probability density will decay like leading to observe AZE instead 
of QZE. As it will be shown below in more detail, note that the effect of introducing N measurements is 
equivalent (regardless of constants) to having , which goes rapidly to zero. 



2.2. Correlation functions and survival probabilities 

Now we shall focus on the quantity central to the discussion in this work: quantum correlation function 
at two different times. Thus, let us consider l^'f^) generically denotes the state of a quantum system at a 
time ti. The unitary time-evolution of this state from ti to ^2 (with t2 > ii) is accounted for the formal 
solution of the time-dependent Schrodinger equation 

I*t2> = U{t2M)\^t,) = e-^^(*-*i)/''|vI/,J. (9) 
The quantum time-correlation function is defined as 

C(t2,ii) = (*tj*t2), (10) 

measuring the correlation existing between the system states at t2 and ti, or, equivalently, after a time 
t — t2 — ti has elapsed (since ti). This second notion also allows us to rewrite (ITUl) as the correlation 
function between the wave function at a time t — t2 — ti and the initial wave function (t = 0), as it follows 
from 

C{t2M) - (*tj*t2> - {■^o\U^{ti)U{t2)\^o) - {^o\U{t2 - ii)|*o) = (*o|*t) = C{t). (11) 
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Another related quantity of interest here is the survival probability, 



P(i2,il)= |(VI/,JVI/, 



Pit). 



(12) 



This quantity indicates how much of the wave function at ti still survives at t2 (in both norm and phase) 
or, equivalently, how much of the initial wave function (also, in norm and phase) survives at a later time 
t = t2 — ti. From now on, concerning the Zeno scenario, we are going to work assuming the second approach, 
although it can be shown that both are equivalent (see Appendix . Taking this into account together with 
the general solution the short-time behavior of P{t) can be readily found. 



P{t)^\{^o\[l- 



iHt 

IT 



2h? 



= 1 - 



fi2 



(13) 



after considering a series expansion up to the second order in t as well as the normalization of V^o- 

In our case, in particular, the fact that the wave packet spreads along time indicates that the quantum 
system becomes more delocalized, this making the corresponding correlation function to decay. This can be 
formally seen by computing the correlation function associated with ([2]), which reads as 



C{t) 



^-Eot^ /2mao{(yo+at)-iEt/h 



-, -1/4 



(14) 



with 



1 



Er^t 1 



1 



- (tan) 



(15) 



(t/2r)2 / 

The exponential in (jl4p only depends on the initial momentum associated with the wave packet centroid, 
but not on its initial position. This is a key point, for the loss of correlation in a wave function displaying a 
translation faster than its spreading rate will mainly arise from the lack of spacial overlapping between its 
values at ti and ^2 (or, equivalently, at and t), rather than to the distortion of its shape (and accumulation 
of phase). However, a relatively slow translational motion will imply that the loss of correlation is mainly 
due to the wave-packet spreading. In this regard, note how the spreading acts as a sort of intrinsic instability, 
which is not related at all with the action of an external potential or a coupling to a surrounding environment, 
but that only comes from the fact that the state describing the system is not stationary (i.e., an energy 
eigenstate of the Hamiltonian, as it would be the case of a plane wave). 

Two scenarios can be thus envisaged to elucidate the mechanisms leading to the natural loss of correlation 
in a quantum system (at this stage, no measurement is assumed). First, consider po — 0, i.e., the wave 
packet only spreads with time, first quadratically and then linearly after the boosting phase [2J], as seen in 
Section 12.11 In this case the correlation function (jl4D reads as 



C{t)=[l + {t/2rf 



-1/4 



(16) 



with 



1 / ^-^ f i 

ipt = (tan) — 

^ 2 ^ l2r 



and the corresponding survival probability (1121) as 

P{t) 



For short times (case (ii) above), ([T^ becomes 



(17) 



(18) 



P{t) « 1 



(19) 



The functional form displayed by (jl9p is the typical quadratic-like decay expected for any general quantum 
state, as it can easily be seen by substituting into the right-hand side of the second equality of with 
Po = 0. Conversely, at very long times. 

Pit) « y , (20) 

i.e., the survival probability decreases monotonically as t~^, in correspondence with the result found above 
for the asymptotic behavior of the probability density. As time evolves, the global phase of the correlation 
function goes from a linear dependence with time (ift ~ —t/Ar) to an asymptotic constant value, Lpoo = — 7r/4. 
Its value thus remains bound at any time between and Lpoo- 

In the more general case of nonzero translational motion for the wave packet (po 7^ 0), the survival 
probability is given by 

pU) = ^ ^-Eot'^ /2mal[l+{t/2Tf] ,21) 

In the short-time limit, this expression reads as 

~ (^1 - ^) ^-Eot^/2v^al_ (22) 

which remarkably stresses the two aforementioned mechanisms competing for the loss of the system corre- 
lation: the spreading of the wave packet and its translational motion. This means that if the translational 
motion is faster than the spreading rate, the wave functions at to and t will not overlap, and P{t) will vanish 
very fast. On the contrary, if the translational motion is relatively slow, the overlapping will be relevant and 
the decay of P{i) will go quadratically with time. In order to express the relationship between spreading 
and translation more explicitly, ((22|) can be expressed in terms of po and p^, i-e.. 



8t2 ; 

Thus, if the translational and spreading motions are such that 



« % (24) 



Vs t 

(actually, it is enough that po/ps ^ l/\/2, since I9>t'^ is already relatively small), then 

2 



P{t) « 1 



1 + 2 



t 



(25) 



which again decays quadratically with time. Otherwise, the decrease of P{t) will be too fast to observe 
either QZE or AZE (see below). Regarding the long-time regime, we find 

P(t) « ^ e-2^=^'/"^o = y e-{po/P.)\ (26) 

which displays the same decay law {t~^) as in the case po = 0, since the argument of the exponential 
function becomes constant. Regarding the phase 5i, it should be mentioned that at short times it depends 
linearly with time, increasing or decreasing depending on which mechanism (translation or spreading) is 
stronger. However, at longer times it approaches asymptotically (also like t"^) the value ipoo regardless of 
which mechanism is the dominant one. 
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3. Quantum Zeno effect and projection operations 



In the standard QZE scenario, a series of von Neumann measurements are performed on tfie system at 
regular intervals of time At. Between two any consecutive measurements the system follows a unitary time- 
evolution according to while each time a measurement takes place (at times t = nAt, with n = 1, 2, . . .) 
the unitarity of the process breaks down and the system quantum state "collapses" into one of the pointer 
states of the measuring device. With this scheme in mind, consider the pointer states are equal to the system 
initial state — in the case we are analyzing here, this type of measurements could consist, for example, of 
a series of diffractions produced by slits with similar transmission properties to the one that generated the 
initial wave function |20]. Thus, after the first measurement the system state will be 

|**=At) = |*o)(*o|*At>, (27) 

which coincides with the initial state, although its amplitude is decreased by a factor (^01^* At)- Each new 
measurement will therefore add a multiplying factor |(^o|^At)P hi the survival probability, which implies 
that it will read as 

^nW= [^i"tT"|(*0|*t-„At)P (28) 

after n measurements, where = |(^o|^At)P- For At sufficiently small, "P^j acquires the form of p^ . 

from which another characteristic time arises, namely the Zeno time [3], defined as 

rz ^ ^. (30) 

In Section 12.11 different stages in the natural evolution of the quantum system were distinguished given 
the ratio between t and the time scale t. The new time scale provided by tz also allows us to distinguish 
between two types of dynamical behavior. For measurements performed at intervals such that <C r^, (^5]) 
holds and the decay of the perturbed correlation function will be relatively slow with respect to the total 
time the system is monitored. Traditionally, this defines the Zeno regime, where the decay of the correlation 
function is said to be inhibited due to the measurements performed on the system. On the contrary, as At 
becomes closer to tz, (HH) does not hold anymore and the decay of the correlation function becomes faster 
than the unperturbed one for finite t. 

In the case of a free Gaussian wave packet with pq — (the system dynamics is only ruled by the wave 
packet spreading), substituting (g]) into (15(11) the Zeno time can be expressed as 



Tz 



= 2V2t, (31) 



which is nearly three times larger than r. According to the standard scenario, provided that At is smaller 
than Tz, one should observe QZE. However, the characteristic time t also plays a key role: as shown below, 
QZE is observable provided that measurements are performed at time intervals much shorter than the time 
scales ruling the wave-packet linear spreading regime. Otherwise, only AZE will be observed. If now we 
consider the more general case, where the free wave packet has an initial momentum (po 7^ 0), a more 
stringent condition is obtained. According to (^5]) — or, equivalently, substituting (jH) into definition ([50]) — , 
we find 

rz - , , (32) 

Vl + 2(poM)2' 

which implies that, in order to observe QZE, the time intervals At between two consecutive measurements 
have to be even shorter (apart from the fact that the condition po/ps ^ 1/V2 should also be satisfied). This 
condition ensures that the wave function at t still has an important overlap with its value at to. 
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Figure 1: (a) Modulus of the time correlation function, |C(t)|, for the unperturbed system (gray line) and three different cases 
with measurements performed at: Ati = 10** <5t = 1 (black), At2 = 10^ St = 0.1 (red), and Ais = 10^ 5t = 0.01 (blue), with 
St = 10~* being the time-step considered in the simulation, (b) and (c) are enlargements of part (a) for times of the order of 
Tz and r, respectively. In the calculations, m = 0.1, (tq = 0.5 and po = 0, which render r = 0.05 and rz = 0.14 (see text for 
details) . 



With the tools developed so far, let us now have a closer look at the QZE and AZE dynamics. Typically, 
these effects are assumed to be quite the opposite. However, we show they constitute the two limits of 
the aforementioned quantum shuffling process. For simplicity and without loss of generality, instead of 
considering the survival probability, in Fig.[IJa) we have plotted the modulus of the time correlation function, 
|C(i)|, against time to monitor the natural (unperturbed) evolution of the wave packet (gray curve) and 
three cases where measurements have been performed at different time intervals At. These intervals have 
been chosen proportional to the time-step 6t (= 10~^ time units) used in the numerical simulation: Ati = 
10" 5t = 1 (black), At2 = 10^ St = 0.1 (red), and Atg = W 8t = 0.01 (blue). Regarding other parameters, 
we have used m = 0.1, cto — 0.5 and po = 0, which make t = 0.05 and tz ~ 0.14. The three color curves 
displayed in Fig. [IJa), which show the action of a set of measurements on the quantum system, behave in 
a similar fashion: they are piecewise functions, each piece being identical to the corresponding one between 

i = and t = At, i.e., to C^^^ = \/^. These curves allow us to illustrate the quantum shuffling process in 
three time regimes which depend on the relationship between r, tz and At: 

(a) For T < Tz < At, the correlation function (see the black curve in Fig. [IJa)) is clearly out of the 
quadratic-like time domain, is convex and therefore the perturbed correlation function always goes 
to zero much faster than the natural decay law (gray curve). This is what we call pure AZE, for the 
correlation function is always decaying below the unperturbed function. 

(b) For T < At < Tz, according to the literature one should observe QZE. However, this is not exactly the 
case. Between r and tz, the correlation function (jl8p displays an inflection point at Tinflx = V^t w 0.071, 
changing from convex to concave. Thus, for At between rinfix and tz, pure AZE is still found due to the 
convexity of the time correlation function. Now, if r < At < Tinfix, the initial falloff of the perturbed 
correlation function is slower, becomes concave and the overall decay gets slower than that associated 
with the unperturbed correlation function (see red curve in Fig.[TJa)). This lasts out for some time, after 
which the perturbed correlation function falls below the unperturbed one (see red curve in Fig.jljb)). It 
is worth stressing here how the decay is indeed faster than in the case of pure AZE, with the quantum 
shuffling making the perturbed correlation function to acquire a seemingly exponential-like shape. 

(c) For At < T, the wave packet is well inside the region where the wave function decay is quadratic-like 
(and concave) and therefore the quantum shuffling produces decays much slower than those observed in 
the unperturbed correlation function as At decreases (see blue curve in Fig. [TJc)). This is commonly 
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Figure 2: (a) Same as Fig. ^a), but showing the envelope II36II superimposed to the corresponding perturbed decay functions: 
Ati = 1 (black), Af2 = 0.1 (red) and Ata = 0.01 (blue). In the figure, the different types of line denote the modulus of 
the time correlation function, |C(t)|, obtained from: the simulation (dotted), the theoretical estimation I I36II (solid), and the 
fitting to a pure exponential function (dashed); to compare with, the unperturbed correlation function is also displayed with 
gray line in panel (a). The decay rates arising from the theoretical estimation are 7J = 25 (black), -y^ = 2.5 (red) and 
■y'^ = 0.25 (blue), while those obtained from the fitting are = 1.224 (black), 7^ = 1.734 (red) and 7^ = 0.249 

(blue), (b) Enlargement of part (a) in the time interval between t = and t = 1. (c) Plot of the difference A(t) between the 
estimated envelope, CAt(i)i a-nd the fitted envelope, in part (a). 



known as QZE. Now, this inhibition of the decay is only apparent; if one considers longer times (see 
blue curve in Fig.[T]^a)), the correlation function is essentially a decreasing exponential, which eventually 
leads the (perturbed) system to decay to zero earlier than its unperturbed counterpart. As it will be 
shown below, these exponential decays can be justified in terms of a sort of Markovianity induced by 
the shuffling process on the system evolution. 

In order to better understand the subtleties behind the quantum shuffling dynamics (and therefore the 



QZE and the AZE), let us focus only on the overall prefactor that appears in 
regime can be written as 

(Ai)2i" 



At 



7(0) 
At 



1 



which in the short-time 



(33) 



This is a discrete function of n, the number of measurements performed up to t„ = nAt, the time at which 
the n-th measurement is carried out. In the limit n — >■ cx), p3p becomes 



-,(00) 

At 



(tn) 



-7Att7i 



(34) 



with the decay rate being 



7Af 



At 



(35) 



as also noted in [2j|. This rate defines another characteristic time, ta* = 7^* ; associated with the falloff of 
the continuous form of 



VAt{t) 



(36) 



(note that this function passes through all the points t„ upon which (l34l) is evaluated). In Fig. [2] we show 
a comparative analysis between the correlation functions |C(t)| displayed in Fig. [T] and their respective 
envelopes, given by CAt{t) = y/VAt{t)', the former are denoted with dotted line and the latter with solid 
line of the same color (again, the gray solid line represents the unperturbed correlation function). The 
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values for the estimated decay rates, given by = for the curves represented, are: = 25 (black), 

72 = 2.5 (red) and 73 = 0.25 (blue). As it can be seen, the agreement between the correlation function 
and its envelope C/^t{t) becomes better as At decreases (see Figs. and (b)), which is in virtue of the 
approximation considered in ((33|) — as Ai increases the behavior of the envelope ((36)) will diverge more 
remarkably with respect to the trend displayed by Pn{t), whereas both will converge as At becomes smaller. 
Thus, while for long intervals At between consecutive measurements the envelope deviates importantly 
from the associated correlation function (see black dotted and solid lines), as At becomes smaller the 
difference between both curves reduces importantly and the relaxation takes longer times (of the order of 
TAt). Nevertheless, for larger values of At one can still perform a fitting of the correlation function to a 
decaying exponential function, C^gj(<) — e^^''**, which renders a qualitatively good overall agreement, as 
can be seen from the corresponding dashed lines in Fig. HJa). Note that the decay rates obtained from 
this fitting are closer to the falloff observed for the corresponding correlations functions (7^ = 1.224, 
72 fit = 1.734, and 73 g^. — 0.249), converging to the estimated value 7' as At decreases (73 « 73). 

From the previous discussion, the idea of a series of sequential measurement acting as a shuffling process, 
wiping out any memory of the system past history, arises in a natural way. Hence, as At becomes smaller, 
\C{t)\ becomes closer to its envelope CAt{t) and therefore to an exponential decay law. Conversely, for larger 
values of At, the system keeps memory of its past evolution for relatively longer periods of time (between two 
consecutive measures), this leading to larger discrepancies between the correlation function and (155)) . Taking 
into account this point of view and getting back to ()28l) , one can notice a remarkable resemblance between 
this expression and a Markov chain of independent processes [13] (which is also inferred from (|34)) ): the 
state after one measurement only depends on the state before it, but not on the previous history or sequence 
until this state is reached. That is, between any two consecutive measurements we have a precise knowledge 
of the probability to find the system in a certain time-dependent state, while, after a measurement, we loose 
any memory on that. Thus, as At becomes smaller, the process becomes fully Markovian, with the time 
correlation function approaching the typical exponential-like decreasing behavior characteristic of this type 
of processes. On the contrary, as At increases, the memory on the past history is kept for longer times, this 
turning the system evolution into non-Markovian, which loses gradually the smooth exponential-like decay 
behavior. Only when a measurement is carried out such memory is suddenly removed, which is the cause 
of the faster (sudden) decays observed in black curve of Fig.[Ha). 

The transition from the non-Markovian to the Markovian regime can be somehow quantified by moni- 
toring along time the distance between the estimated envelope, CAt{t), and the fitted envelope, CAt,fit{t), 

A{t)=CAt{t)^CAtMt), (37) 

which is plotted in Fig.[5fc) for the three cases of At considered. Thus, as At becomes smaller, we approach 
an exponential decay law and A(t) goes to zero for any time (see blue curve in the figure), this being 
the signature of Markovianity. On the contrary, if the time evolution is not Markovian, as time increases 
and the system keeps memory for longer times, the value of A(t) displays important deviations from zero 
(see black and red curves in the figure). These deviations mainly concentrate on the short and medium 
term dynamics, where values are relatively large to be remarkable. Nonetheless, analogously, one could 
also display the relative ratio between the two correlation functions, which would indicate or not the trend 
toward Markovianity in the long-time (asymptotic) regime. 

4. Conclusions 

By assuming that the QZE inhibits the evolution of an unstable quantum system, one might also be 
tempted to think that its coherence is also preserved, while the AZE would lead to the opposite effect in its 
way through faster system decays. In order to better understand these effects, here we have focused directly 
on the bare system, i.e., no external potentials or surrounding environments acting on the quantum system 
have been assumed. This has allowed us to elucidate the conditions under which such effects take place in 
relation to the intrinsic time-scales characterizing the isolated system, which have been shown to play an 
important role. Furthermore, by means of this analysis, we have also shown that both QZE and AZE are 
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indeed two instances of a more general effect, namely a quantum shuffling process, which eventually leads the 
system to display a Markovian-like evolution and its correlation function to follow an exponential decay law 
as the interval between measurements decreases. Within this scenario, the QZE dynamics can be regarded 
as a regime where any trace of knowledge on the initial system state is lost due to a rapid shuffling, while 
in the long-time regime the correlation function would fall to zero faster than the unperturbed one. The 
apparent contradiction with the traditional no-evolution scenario can be explained very easily: Since one 
often cares only about the short-time dynamics, the long-time dynamics is usually completely neglected. 
In other words, the time during which the system dynamics is usually studied is relatively small compared 
to the Markovian time-scale induced by the continuous measurement process and, therefore, one assumes 
nearly stationary dynamics. 
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A. Alternative Zeno scenario 

In the Zeno scenario considered above, establishing a direct analogy with Zeno's arrow, the wave packet 
plays the role of a quantum arrow, but with the particularity that this arrow slows down until its evolution 
is frozen by means of a series of measurements. This is the scenario traditionally considered Q. However, 
a more direct analogy with Zeno's arrow paradox can be established if it is assumed that the wave packet 
is always in motion and the measurements are just like photographs indicating the particular instant from 



which the correlation has to be computed 18|. Because of the actions on the wave packet in relation to 
what a measurement is considered in each, we can call these two situations as: 

(a) The stopping- arrow scenario, where the wave packet is "collapsed" or "stopped" after each measurement. 

(b) The steady-arrow scenario, where the wave packet time-evolution never stops, but the computation of 
the correlation function is reset after each (photograph-like) measurement — like in a stop-motion or 
stop-action movie. 

It can be shown that both scenarios are equivalent with the exception of a lost time-dependent phase 
in the latter. In order to prove this statement, let us start by considering the wave function is now left to 
freely evolve in time. Following the idea behind this scenario, the survival probability is monitored in time 
by computing the overlapping of the wave function at a time t with its value at successive times to, ti, t2, 
with tn = nAt, i.e., 

(^t„-J*t>, (38) 
where tn-i < t < tn- Note here the direct analogy with Zeno's arrow, where at each instant we are observing 
the arrow steady at a different space position, but without freezing its motion. Taking this into account, for 
< t < ti, we have 

P(t) = |(vI/o|**)p. (39) 

Now, if ti < < < t2, 

P(t)=«i|(vl/,jvl/,)|2, (40) 
while the wave function for the same interval will be 

|*t) = V^\-^t), (41) 

i.e., the evolved wave function, but with a prefactor which ensures the matching of the different of P{t) 
before and after t — ti. Following the same argumentation, after the n — 1 measurement, 

P{t) = {KZlo^k)\{-^t^^A^t)\\ (42) 
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with tn-l < t < tn- 

As it can be seen, by means of this procedure the wave function is never altered (which always keeps 
evolving according to the Schrodinger equation), but only its relative amplitude. So, the key issue here is 
the attenuation factor a, which can be evaluated as follows. For < t < t„, we note that 

{^'t„_A'^t) = («'o|e*-^("-i)^*/'"'e~'-^*/'"'|^'o) = (*o|e'-^'*^'""^''^*'^''|*o). (43) 

Therefore, the attenuation factor for the interval i„_i < t < t„ should come from the overlapping of the 
wave function at the times when the two previous measurements were performed, i.e., 

a„_i = K^t^jM/t^Jp |(*o|e'^^*/''|vE'o)p. (44) 

This expression can be Taylor expanded to second order in At (assuming At is small enough) taking into 
account that 

2 

(*o|*At) ~ 1 - ^ (*o|i?|*o) + ^ (^ol-ff'l^'o), (45) 
where we have assumed the wave function is initially normalized. This renders 

a„ = K^ol^-t,)!' « 1 - ^ {{^o\H'\^o) - (*o|i?|*o)') = 1 - ^ (Ai^)^ (46) 

which is valid for any n, since it does not depend explicitly on the particular time at which the measurement 
is made. Therefore, (H^ can be expressed as 

Pit)^ar'm„-A^t)\', (47) 

which is formally equivalent to (|28p after performing n measurements, although it describes a completely 
different physics [l8[. 
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